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4.4 Buoyant plume from a steady heat source 

[Reference] : 

Gebhart, et. al. (Jalluria, Maharjan, Saammakia), 
Buoyancy-induced Flows and Transport, 1988, 
Hemisphere Publishing Corporation 

Let T = T — Too = temperature variation where Tiufty is a constant (no ambient 
stratification). For a strong enough heat source, we expect boundary layer behavior, 
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The boundary layer equations are 
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The centerline r = is an axis of symmetry, 

du df ^ 

Far outside the plume r — )■ oo 

w ^ and T ^ Too, (f ^ 0) (4.4.5) 

Rewrite (4.4.3) as 
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after using continuity. Now integrating the last equation from r = to r = oo 
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Therefore 
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Using the boundary conditions, we get or 
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2'KruTdr — constant 



Note that 



/ 2Tirdr upCT = rate of buoyancy flux 
= rate of heat flux 

= (5 (given rate of heat release at a; = 0) 



Therefore, 



(5= / 2'KrdrfmCT 
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This is a boundary condition. 

Let the stream function be defined by 
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dip 



dr ' dx 
(4.4.1) is automatically satisfied. From the momentum equation: 
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From the energy equation 
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and from the buoyancy flux condition 
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Try a similarity solution with the one-parameter transformation 

x-A"x*, r = AV, V = A'V*, f = X^T* 

Prom (4.4.10), 

^2c-ih-a ^ ^2c-46-a ^ ^ (4.4.13) 

from (4.4.11) 

^c+d-2b-a ^ (4.4.14) 

and from (4.4.12) 

A'^+'' = 1 (4.4.15) 

Prom these three equations we get 

c b_l d_ 

a a Z a 

Th similarity variable can be taken to be 

V^^j, (4.4.16) 



^1/2 

and the similarity solution to be 

'ilj^xF{'q) (4.4.17) 

and 

f^x-^G{ri) (4.4.18) 

After a lot of algebra, and noting 

drj 1 drj 1 r Irl rj 

dr ' dx 2 x^/"^ 2 x^/^ x 2x 

we get from (4.4.10) 

uF'" + (^^^ (F-u)+ gprjG = (4.4.19) 

and from (4.4.11) 

k{r]G'y + {FCy = (4.4.20) 

Before integrating, let us normalize : 

rj^afj, F^jF, G^aG. (4.4.21) 

It follows from (4.4.19) that 

^,F"' + ^ {—\ (7F - I/) + gpaafiG = (4.4.22) 
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where prime denotes d/dr]. Setting ^ = v and 

3 = g(3aa 



which relates a and a, 



a = -— (4.4.23) 
we get 

i^'" + {F-l)+f]G = (4.4.24) 

Similar normalization of (4.4.20) gives 

^(77(5')'+^(F(5)' = (4.4.25) 



. which can be simplified to 
. where 



{fjCy + Pr{FG)' = (4.4.26) 



1/ 



P^ = - = Prandtl Number (4.4.27) 

rC 

For water u = lO^^cm^/s, /c = 1.42cm^/s, hence Pr = 7. For air i/ = 0.145cm^/s, /c = 
0.202cmVs, hence Pr = 0.75. 

We now integrate (4.4.26)to give 

fjG' + PrFG = constant 

Since ip{x, 0) = 0, we must have P(0) = ; the constant above is zero. 

fiG' + PrFG = (4.4.28) 

Equation (4.4.28) can be written 

G' „ P d In G „ P 
= -Pr-, or = -Pr- 

G rj drj rj 



constant 



In G = -Pr / -d-q + ' 
Jo 77 

G{f}) = G'(O) exp l^-P^ £ ^dn^ (4.4.29) 

Substituting Eqn. (4.4.29) into Eqn. (4.4.24), the resulting equation for P must be integrated 
numerically. 

Now let us find the boundary condtions for P or P. 
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Eqn. (4.4.8 ) becomes 



zTrpC 7o \r or J X Jo r x Jo Jo 

(4.4.30) 

Therefore, 

Q_ 
2Tr pCua 

Let us choose 

1 (4.4.32) 



/ dfjF'G=-^^ (4.4.31) 

Q 

2t{ pCva 
so that 

/■oo _ _ 

/ dfjF'G=l (4.4.33) 

JO 

is the boundary condition for F and G. Now (4.4.32) defines a, the scale of G. Note that 
larger Q implies larger a and smaller a. Thus a stronger heat source leads to a greater 
center line temperature and a thinner plume. Also, 

It — >■ as r — )■ oo 

hence 

1 , F' u F' 

U = -iPr = = >-0, as T] ^ T] ^ OO 

r r] r] 



The radial velocity is, in general 



1 , 1 F'^ 



Since 

we must have. 
Clearly 



V ^ as ?7 ^ 0, 
F(0) = 0. 
F' 

— 0, and F{f}) = as fy (4.4.34) 

The numerical results by Mollendorf & Gelbart, 1974, are shown in Figs. 4.4.1, for 
various Prandtl numbers. A schlierian photograph due to Gebhart (copied from Van Dyke 
An Album of Fluid Motion) is hown in Figure fig: plume VD. 

Remark : 

Idi^ F' ( X ^, I 
V 



r dr ?7 V r x^/^ 
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Along the centerline u{x, 0) = ("^j^ = constant depending on P^- Why? Buoyancy acceler- 
ation is counteracted by entrainment. 

Remark : Let the radius of the plume be a which varies as 

a ~ x'/' 

This is consistent with the behavior that u x^, and T ~ since 

a'^uf = Q 

On the other hand the mass flux rate is 

ua^ ~ X 

and the momentum flux rate is 

2 2 

u a ^ X 

hence both approach zero at the source. Thus a plume is the result of energy source, not of 
mass or momentum. 




(Adaptsd fronj MoiSendorf and Getthatt. J 974}. 



